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ABSTRACT 

It is now accepted that accretion onto classical T Tauri stars is controlled by the 
stellar magnetosphere, yet to date most accretion models have assumed that their 
magnetic fields are dipolar. By considering a simple steady state accretion model 
with both dipolar and complex magnetic fields we find a correlation between mass 
accretion rate and stellar mass of the form M oc M", with our results consistent 
within observed scatter. For any particular stellar mass there can be several orders of 
magnitude difference in the mass accretion rate, with accretion filling factors of a few 
percent. We demonstrate that the field geometry has a significant effect in controlling 
the location and distribution of hot spots, formed on the stellar surface from the high 
velocity impact of accreting material. We find that hot spots are often at mid to low 
latitudes, in contrast to what is expected for accretion to dipolar fields, and that 
particularly for higher mass stars, the accretion flow is predominantly carried by open 
field lines. 

Key words: Stars: pre-main sequence - Stars: magnetic fields - Stars: spots - Stars: 
formation - Stars: coronae - Stars: low mass, brown dwarfs 



; 1 INTRODUCTION 

. Classical T Tauri stars (CTTSs) are young, low mass, pre- 
' main sequence stars that are actively accreting from a sur- 
. rounding disc whi c h is t he eventual birth-place of planets. 
' lUchida fc Shibatal lll984f) suggested that the magnetic field 
' of a CTTS disrupts the inner disc. In the early 1990s several 
magnetosphcric a ccretion models w ere developed (iKonigll 
[1991; Collier Ca meron fc Campbelll [l993: Shu ct al.l 1199411 
where material is lifted from the disc plane and is channelled 
along dipolar magnetic field lines onto the star, terminating 
in a shock at the photosphere. In an idealised model of a 
CTTS's magnetic field there are closed field lines close to 
the star that contain the X-ray emitting corona, whilst at 
larger radii, there are closed field lines which thread the cir- 
cumstellar disc. It is along this latter set of field lines that 
accretion may proceed. There are also regions of open field 
which carry outfiows in the form of a wind, and in some 
cases, as large coUimated bipolar jets. 

Magnetospheric accretion models assume that CTTSs 
possess magnetic fields that are strong enough to dis- 
rupt the disc at a distance of a few stellar radii. Such 
strong fields have been detected in a number of systems 
using a variety of techniques. Average surface fields of 
l-3kG have been detected most successfully by exploit- 
ing the Zeeman effect, both through Zeeman broaden- 
ing (e.g. Ijohns-Krull. Valenti fc Koreskd ll999lJl and from 
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the circular polarisation of lines whic h are sensitive to 
the p r esence of a magnetic field (e.g. |johnsrKrulLet_al 
|l999al : ISvmington et al.ll2005l: lOaou. Johns-KruU fc Valenti 
I2OO6I) . Field detections have also been made from the 
increase in line equivalent width iBasri. Marcv fc Valentil 
ll992l:lGuent~ et al .ll999^ and also from electron cyclotron 
maser emission, a coherent emission process from mildly rel- 
ativistic electrons t rapped inside fiux tubes close to the star 
JSmith et alj|2003ll . The mean magnetic field strengths de- 
tected so far appear to be ro ughly constant across all stars 
iValenti fc Johns- Krullll2004ri . 

Traditionally magnetospheric accretion models have as- 
sumed the CTTSs have dipolar magnetic fields. Dipole fields 
(or inclined dipole fields) have been successively used to ex- 
plain some of the observations of C TTSs (e.g. the photopo - 
larimetric variability of A A Tau. lO 'Sullivan et alJ l2005| l. 
but fail to account for others. fValenti <^^oIm^<nfi]|l20M) 
present magnetic field measurements for a number of stars, 
and despite detecting strong average surface fields from Zee- 
man broadening, often measurements of the longitudinal 
(line-of-sight) field component (obtained from photospheric 
lines) are consistent with no net circular polarisation. This 
can be interpreted as there being many regions of opposite 
polarity on the stellar surface, giving rise to oppositely po- 
larised signals which cancel each other out giving a net polar- 
isation signal of zero. This suggests that CTTSs have mag- 
netic fields which are highly complex, particularly c l ose to 
the stellar surface; however, as lValenti fc Johns-Krull (|20^ 
point out, as the higher order multi-pole field components 
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will drop off quickly with distance from the star, the dipole 
component may still remain dominant at the inner edge of 
the disc. Also their measurements of the circular polarisa- 
tion of the Hel 5876A emission line (believed to form in the 
base of accretion columns) are well fitted by a simple model 
of a single magnetic spot on the surface of the star, suggest- 
ing that the accreting field may be well ordered, despite the 
surface field being complex. 

The fractional surface area of a CTTS which is 
covered in hot spots, the accretion filling factor face, 
is inferred from obs ervations to be small; typically of 
order one percent ij ^izeroll^t ^ l2003l: ICalvet et alJ 
| 2004l: IValenti fc Johns-KrulJ l2004l : ISvmington et al.l 120051 
iMu^eroU^^^n 120051? ! Dipolar magnetic field mod- 
els predict accretion filling factors which are too 
large. This, c ombined with the polarisation results, led 
IJohns-Krullfc Ga fford ( 200i) to generalise the Shu X-wind 
model JShu et alJll994h to include multipolar, rather than 
dipolar, magnetic fields. With the assumption that the av- 
erage surface field strength does not vary much from star 
to star the generalised Shu X-wind model predicts a corre- 
lation between the stellar and accretion parameters of the 
form Rtfacc oc {M^M ProtY^'^ , a prediction that matches 
observations reasonably well. 

In this paper we present a model of the accretion pro- 
cess using both dipolar and complex magnetic fields. We 
apply our model to a large sample of pre-main sequence 
stars ob tained from the Cha ndra Orion Ultradeep Project 
(COUP; iGetman et aIll2005^ . in order to test if our model 
can reproduce the observed correlation between mass accre- 
tion rate and stellar mass. An in c rease in M with M« was 
ori ginally noted b y iRe buU et a l] l|2000 () and su bsequently 
bv i White fc Ghed i200lll and lRebuU et alj ll2002ll . The cor- 
relation was then found to extend to very lo w mass ob- 
jects and accreting brown dwarfs bv lWhite fc Basri (,2003) 
and iMuzeroUe et alJ i2003l) , and to the higher mass, inter- 
mediate mass T Tauri stars, by ICalvet et al.l (|200^. Fur- 
ther lo w mass data has recently been adde d bv | Natta et alJ 
liooi), 'Mohantv. Jav awardhana fc Basril ^2005^ and bv 
rMuzcroUc ct a l. ( 2005) who obtain a correlation of the form 
M oc M1'^ , with as much as three orders of magnitude 
scatter in the measured mass accretion rat e at an y par- 
ticular stellar mass. However, ICalvet et alJ (l2004tl point 
out that due to a bias against the detection of higher 
mass stars with lower mass accretion rates, the power may 
be less than 2.1. Further data for accreting stars in the 
p— Ophiuchus star forming region has recently been added 
bv lNatta. Testi fc RandichI (l2006ll . 

The physical origin of the correlation between M* and 
M, and the large scatter in measured M values, is not 
clear; however several ideas have been put forward. First, 
increased X-r ay emission in higher mass T Tauri stars 
jPreibisch et a l. 2005: Jardine ct al. 20063 may cause an in- 
crease in disc ionisation, leading to a more efficient magne- 
torot ational instability a nd therefo re a higher mas s accre tion 
rate JCalvet et alj|2004ll . Second. IPadoan et alj (|20oi) ar- 
gue that the correlation M oc arises from Bondi-Hoyle 
accretion, with the star-disc system gathering mass as it 
moves through the parent cloud. In their model the observed 
scatter in M arises from var iations in stella r veloc ities, gas 
densities and sound speeds. iMohantv et alJ i2005f) provide 
a detailed discussion of both of these suggestions. Third, 
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[Alexander fc Armitagd i200(i) suggest that the correlation 
may arise from variations in the disc initial conditions com- 
bined with the resulting viscous evolution of the disc. In 
their model they assume that the initial disc mass scales 
linearly with the stellar mass, Md oc M«, which, upon mak- 
ing this assumption, eventually leads them to the conclusion 
that brown dwarfs (the lowest mass accretors) should have 
discs which are larger than higher mass accretors. However, 
if it is the case that the initial disc mass increases more 
steeply with stellar mass, Md oc , then the stellar mass 
- accretion rate correlation can be reproduced with smaller 
brown d warf discs of low mass (o f order one Jupiter mass). 
Thus the [Alexander fc Armitagd i2006l) suggestion, if cor- 
rect, will soon be d irectly verifiable by observations. Fourth, 

iNatta et alj (1200(1) suggest that the large scatter in the cor- 
relation between M and Af« may arise from the infiuence 
of close companion stars, or by time variable accreti on. It 
should however be noted that IClarke fc Pringl^ (l2006f) take 
a more conservative view by demonstrating that a steep cor- 
relation between M and M« may arise as a consequence of 
detection/selection limitations, and as such M oc is per- 
haps not a true representation of the correlation between M 
and M* . 

In §2 we describe how magnetic fields are extrapolated 
from observed surface magnetograms. In §3 we consider ac- 
cretion onto an aligned, and then a tilted dipole field, to 
develop a simple steady state accretion model and to inves- 
tigate how tilting the field affects the mass accretion rate. In 
§4 these ideas are extended by considering magnetic fields 
with a realistic degree of complexity and we apply our accre- 
tion model to study the correlation between mass accretion 
rate and stellar mass, whilst §5 contains our conclusions. 



2 REALISTIC MAGNETIC FIELDS 

From Zeeman-Doppler images it is possible to extrapolate 
stellar magnetic fields by assuming that the field is po- 
tential. At the moment we do not have the necessary ob- 
servations of CTTSs, but we do have for the solar like 
stars LQ Hya and AB Dor ([D onati fc Collier Cameroiil 
| l997t iDonati et alJll997l : iDonatil 119991 : iDonati et alJll999l : 
bonati et al.ll2003f). which have different fi eld topologies 
(Ijardinc. Collier Cameron fc Donatil l2002al: iHussain et alJ 
12002.: .Mclvor et alj|2003l : iMcIvor et al J 120041) . Using their 
field structures as an example we can adjust the stellar pa- 
rameters (mass, radius and rotation period) to construct a 
simple model of a CTTS, surrounded by a thin accretion 
disc. 

The method for extrapola t ing m agnetic fields follows 
that employed bv l Jardine et alJ(l2002a^ . Assuming the mag- 
netic field B is potential, or current-free, then V x B = 0. 
This condition is satisfied by writing the field in terms of a 
scalar flux function ^, such that B = — V^'. Thus in order 
to ensure that the field is divergence- free (V • B = 0), ^I" 
must satisfy Laplace's equation, — 0; the solution of 

which is a linear combination of spherical harmonics, 

JV I 

* = ^ ^ [aimv + bi^r-^'+''>] P^™(e)e""^ (1) 

l — l m — — l 

where Pim denote the associated Legendre functions. It then 
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follows that the magnetic field components at any point 
(r, e, (j)) are, 



Table 1. Data for CTTSs from lValenti fc Johns-KrulJ Eog 
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N I 



-B^ = -2^ 2_^\aimr +bijnr '\ — —-ime ^. (4) 



The coefficients ai^ and bim are determined from the radial 
field at the stellar surface obtained from Zeeman-Doppler 
maps and also by assuming that at some height Rs above 
the surface (known as the source surface) the field be- 
comes radial and hence Bg{Rs) = 0, emulating the effect 
of the corona blowing open fi eld lines to form a stellar wind 
jAltschuler fc Newkirk^ll969^ . In order to extrapolate the 
field we used a modified version of a code originally devel- 
oped bv lvan Balleeooiien. Cartledee fc Priesti 1179981) . 



2.1 Coronal extent 

We determine the maximum possible extent of the corona 
(which is the extent of the source surface) by determining 
the maodmum radius at which a magnetic field could contain 
the coronal gas. Since a dipole field falls off with radius most 
slowly, we use this to set the source surface. For a given sur- 
face magnetogram we calculate the dipole field that has the 
same average field strength. We then need to calculate the 
hydrostatic pressure along each field line. For an isothermal 
corona and assuming that the plasma along the field is in 
hydrostatic equilibrium then. 



Ps 



■■ po exp , -2 



(5) 



where Cs is the isothermal sound speed and Qs the component 
of the effective gravity along the field line such that, Ps = g- 
B/|B|, Po is the gas pressure at a field line foot point and ps 
the pressure at some point along the field line. The effective 
gravity in spherical coordinates for a star with rotation rate 
u} is, 



GM« 2 . 2 
— \- iL) r sm ( 
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We can then calculate how the plasma j3, the ratio of gas 
to magnetic pressure, changes along each field line. If at 
any point along a field line (3 > 1 then we assume that the 
field line is blown open. This effect is incorporated into our 
model by setting the coronal (gas) pressure to zero whenever 
it exceeds the magnetic pressure (/3 > 1). We also set the 
coronal pressure to zero for open field lines, which have one 
foot point on the star and one at infinity. The gas pressure, 
and therefore the plasma /3, is dependent up on the choice 
of Po which is a free parameter of our model. Ijardine et alJ 
i2006l) provide a detailed explanation of how po, the coronal 
base (gas) pressure, can be scaled to the magnetic pressure 
at a field line foot point, so we provide only an outline here. 
We assume that the base pressure is proportional to the 
magnetic pressure then po = KBq, a technique which has 
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been used successfully to calculate mean coronal densities 
and X-ray emi ssion meas ures for the Sun and other main 
sequence stars ijardine et al.. .2002al lbl). By varying the con- 
stant K we can raise or lower the overall gas pressure along 
field line loops. If the value of K is large many field lines 
would be blown open and the corona would be compact, 
whilst if the value of K is small then the magnetic field 
is able to contain more of the coronal gas. The extent of 
the corona therefore depends both on the value of K and 
also on Bo which is determined directly from surface mag- 
netograms. For an observed surface magnetogram the base 
magnetic pressure Bo varies across the stellar surface, and 
as such so does the base pressure po at field li ne foot points. 
By co nsidering stars from the COUP dataset [jardine et alJ 
i2006i) obtain the value of K which results in the best fit to 
observed X-ray emission measures, for a given surface mag- 
netogram (see their Table 1). We have adopted the same 
values in this paper. We then make a conservative estimate 
of the size of a star's corona by calculating the largest radial 
distance at which a dipole field line would remain closed, 
which we refer to as the source surface Ra. 



2.2 Coronal stripping 

Lower mass stars, which have small surface gravities and 
therefore large pressure scale heights, typically have more 
extended coronae which would naturally extend beyond the 
corotation radius. Closed field lines threading the disc be- 
yond corotation would quickly be wrapped up and sheared 
open. Therefore if the maximum extent of a star's corona 
is greater than the corotation radius then we set it to be 
the corotation radius instead. Higher mass stars, with their 
larger surface gravities, typically have more compact coro- 
nae which may not extend as far as the corotation radius. 
Therefore if we assume that the disc is truncated at corota- 
tion, then accretion proceeds from the inne r edge of the disc 
along radial open field lines. Ijardine et alJ i2006l) provide a 
discussion about the extent of T Tauri coron ae rela t ive to 
corotation radii. It is also worth noting that ISafieil lll998l) 
criticises current magnetospheric accretion models for not 
including the effects of a stellar corona. He argues that the 
inclusion of a realistic corona blows open most of the closed 
field with the eventual net effect being that the disc would 
extend closer to the star. However, in our model there are 
open field lines threading the disc at corotation and it is 
therefore reasonable to assume that they are able to carry 
accretion flows. 



2.3 Field extrapolations 

The initial field extrapolation yields regions of open and 
closed field lines. Of these, some intersect the disc and may 
be actively accreting. For closed field lines which do not in- 
tersect the disc it is possible to calculate the X-ray emission 
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(a) (b) 

Figure 1. Field lines which could support accretion flows for a model of a CTTS with a field topology that resembles (a) LQ Hya, 
obtained using the DF Tau parameters from Table HI and (b) AB Dor using the CY Tau parameters. The stellar surface is coloured to 
show the strength of the radial component of the field, with red representing IkG and black -IkG. Field lines have been drawn from the 
corotation radius. For the lower mass star, DF Tau, the natural extent of its corona would be beyond corotation and therefore there is a 
mixture of open and closed field lines threading the disc at Rco- The higher mass star, CY Tau, has a more compact corona and material 
flows along open field lines from corotation. 



measure in the same way as Ijardine et alJ i2006l) . To de- 
termine if a field line can accrete we find where it threads 
the disc and calculate if the effective gravity along the path 
of the field line points inwards, towards the star. From this 
subset of field lines we select those which have /3 < 1 along 
their length. In other words, for any given solid angle we 
assume that accretion can occur along the first field line 
within the corotation radius which is able to contain the 
coronal plasma. We assume that the loading of disc mate- 
rial onto the field lines is infinitely efficient, such that the 
first field line at any azimuth which satisfies the accretion 
conditions will accrete, and that field lines interior to this 
are shielded from the accretion flow. We also assume that 
the accreting field is static and is therefore not distorted by 
the disc or by the process of accretion. In §3.3 we consider 
in more detail how to determine which field lines are able to 
support accretion flows, in order to calculate mass accretion 
rates and accretion flUing factors. 

Fig. shows the first set of field lines which may be 
accreting, obtained by surrounding the field extrapolations 
of LQ Hya and AB Dor with a thin wedge-shaped accretion 
disc, with an opening angle of approximately 10°. In §3 we 
develop a model for isothermal accretion flows where mate- 
rial leaves the disc at a low subsonic speed, but arrives at 
the star with a large supersonic speed. Not all of the field 
lines in Fig.^are capable of supporting such accretion flows, 
and instead represent the maximum possible set of held lines 
which may be accreting. We assume a coronal temperature 
of lOMK and obtain the gas press ure at the base of eac h field 
line as discussed in §2.1 and bv lJardine et al.) l|200^. The 
natural extent of the corona of DF Tau would be beyond 
the corotation radius and therefore accretion occurs along a 



mixture of closed and open field lines from corotation. One 
suggestion for how accretion may proceed along open field 
lines is that an open field line which stretches out into the 
disc, may reconnect with another open field line for long 
enough for accretion to occur, only to be sheared open once 
again. This is of particular importance for the higher mass 
stars, such as CY Tau, where in some cases we find that 
the inner edge of the disc is sitting in a reservoir of radial 
open field lines. This may have important implications for 
the transfer of torques between the disc and star. However 
more work is needed here in order to develop models for 
accretion along open field lines. 

These field extrapolations suggest that accretion may 
occur along field lines that have very different geometries. 
Indeed, a substantial fraction of the total mass accretion 
rate may be carried on open field lines. Before developing a 
detailed model of the mass accretion process, however, we 
first consider the simple case of a tilted dipole. While this is 
an idealisation of the true stellar field it allows us to clarify 
the role that the geometry of the field may have in governing 
the mass accretion process. 



3 ACCRETION TO A DIPOLE 

We have constructed two simple analytic models as sketched 
in Fig. |5| The first case is for a star with a dipolar field 
with the dipole moment /x aligned with the stellar rotation 
axis f2. In standard spherical coordinates this field may be 
described as, 

B=(^cose,iisine,0), (7) 
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Figure 2. An aligned and tilted dipolc field geometry. The 
aligned dipole (left) with a field line in the star's meridional plane, 
with the dipole moment fj, aligned with the stellar rotation axis f2, 
and the perpendicular dipole (right) with a field line in the star's 
equatorial plane, with /x perpendicular to O. The average s urface 
field strength matches that considered by I.Tardine et al.l J200fih 
with yellow (blue) representing the positive (negative) magnetic 
pole. 



a scenario that allows us to model accretion flows along field 
lines in the star's meridional plane. If we then take this field 
structure and tilt it by 7r/2 radians such that /i now lies in 
the star's equatorial plane, perpendicular to £7, then those 
field lines which ran north-south in the meridional plane, 
now lie east-west in the equatorial plane, with, 



/ 2u u 
B = — cos (h, 0, sm ( 



(8) 



Throughout we shall refer to these cases as the perpendic- 
ular dipole for the tilted dipole field and the aligned dipole 
for the aligned dipole field. To do this we consider steady 
isothermal accretion fiows from a thin accretion disc ori- 
ented such that the disc normal is parallel to the stellar 
rotation axis. An initial sonic Mach number, M, is ascribed 
to the accreting material. We then calculate the pressure 
and velocity profiles, relative to arbitrary initial conditions 
defined at the disc plane. We calculate the ratio of pressure 
p at each point along a field line, relative to that at the disc, 
Pd', and then from this we calculate how the Mach number 
of the flow changes along the field. 

The path of a field line may be described by 



Br 

dr 



Bg_^ 

rd6 r sin 



(9) 



For the perpendicular dipole, where 6 = n/2 for all field lines 
that pass through the disc, it is quickly established that. 



(10) 



where ^' is a constant along a particular field line, such that 
different values of *]/ correspond to different field lines. For 
the perpendicular dipole case, at r — Rd, the maximum 
radial extent of the field line, (j) = 7r/2; thus, ^ = 1/Rd- 
Using this result the magnitude of the magnetic field at a 
point along a field line a distance r from the centre of the 
star, relative to that at the disc, may be obtained from ©, 



B_ 
B'd 



Rd 



(11) 



where Bd — B [r = Rd). An identical expression can be de- 
rived for the aligned dipole case. 

3.1 Steady isothermal accretion flows 

The momentum equation for a steady inviscid fiow along a 
flux tube is 



p(v- V)v = - V p-F 



B2 



(B • V)B-hpg,„-2ptJxv, (12) 



where the symbols have their usual meaning with ge/ / being 
the sum of the gravitational and centrifugal accelerations. In 
a frame of reference rotating with the the star the effective 
gravity is. 



: g — X (tj X r) , 



(13) 



The component of the momentum equation along a field line 
is then, 



dp 



(14) 



since the Coriolis term {—2puj x v) does not contribute for 
flows along the field, and where s(r) is a unit vector along 
the path of the fleld line. Throughout terms with a sub- 
script d will denote quantities defined at the disc; for exam- 
ple pd,Pd,Vd and Bd are respectively the density, pressure, 
velocity along the field and the magnetic field strength as 
defined at the plane of the disc, a radial distance Rd from 
the centre of the star. Integrating equation (I14I I from the 
disc plane to some position along the field line at a distance 
r from the stellar centre, and using the isothermal equation 
of state for an ideal gas, p = pCg, gives. 



In 



2^ 



Seff 



■ sds 



(15) 



(pvA) = 
(BA) = 0, 



where Cs is the isothermal sound speed. If we assume that 
both mass and magnetic flux are conserved along each flux 
tube (of cross-sectional area A), then the fiow must satisfy 

d_ 

ds 
d_ 

ds 

which may be expressed equivalently as, 

PV _ PdVd 

B ^ Bd 

By combining 11511 with (I18II a relation for the pressure struc- 
ture along an accreting field line can be established, 

and also directly from I18II an expression for the velocity 
structure. 



= const, 



(16) 
(17) 

(18) 



M 



Pd B 



(20) 



P Bd' 

where in both cases Ml = I'd/cs denotes the initial sonic 
Mach number at which the accretion flow leaves the plane 
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of the disc. It is then possible to find the pressure at each 
point along a field line, relative to the pressure at the disc 
{p/pd), by finding the roots of Hf 9|l . Once these roots have 
been found the velocity profile can be obtained from 12UII 
by calculating how the Mach number of the flow varies as 
material moves from the disc to the star. 



3.2 Pressure and velocity profiles 

For the perpendicular dipole, in the star's equatorial plane, 
the effective gravity has only a radial component, 

GA'L 



Seff 



+ LU r r. 



(21) 



Taking the component of the effective gravity along the field, 
that is along a path parameterised by § = B/B, substituting 
into II19II , and using lllll gives an expression for the pressure 
structure along equatorial field lines. 



Pd 



$c (Ri - r") = 0, (22) 



where both r and Rd are measured in units of the stellar 
radius J?, and $g and $c are the surface ratios of the grav- 
itational and centrifugal energies to the thermal energies, 

GA'L 



1 /ioR, 



(23) 
(24) 



The roots of i'2'21 give the pressure at some point along a 
field line loop which is a radial distance r from the stellar 
centre. 

For the aligned dipole, in the star's meridional plane, 
the effective gravity has both an r and 6 component, 

GAL 



Seff 



+ uj^r sin^ 9] f + (uj'^r sin 6 cos 6) 6. (25) 



Following an identical argument to that above it can be 
established that for accretion in the star's meridional plane, 
the pressure function (If 91 becomes. 



""Si + 2 



5-'(^)"(- 



3r ^^ f Pd 

Rd 



Rd r 



Rd 



■ 0. 



(26) 



The only difference from the perpendicular dipole is in the 
final term. 

For a CTTS with a mass of 0.5Mq, radius 2Ro and 
a rotation period of 7 days we have calculated the pressure 
and velocity structure along accreting dipole field lines, for a 
range of accretion flow temperatures, starting radii and ini- 
tial sonic Mach numbers. Figs. |3(ei)| and [Sj^b)! show a typical 
pressure and velocity profile for the perpendicular dipole, 
whilst those for the aligned dipole are qualitatively similar. 
The pressure profile shows how the ratio p/pd, where p is the 
pressure along the field line and pd the pressure at the disc, 
varies as the flow moves from the disc to the star (plotted 
logarithmically for clarity). The velocity profile shows how 
the Mach number of the flow changes along the field line. 



For different accretion flow temperatures and starting radii 
the resulting profiles are similar, except in a few select cases, 
as discussed in the next section. Fig. |H] is for an accretion 
flow leaving the disc at Rd — 6. OR*, which is approximately 
the equatorial corotation radius Rco where. 



Rco — I 



1/3 



(27) 



and for an accretion flow temperature of fO^K. This is at 
least one order of magnitude higher than what is believed 
to be typical for accretion in CTTSs, but a higher temper- 
ature has been selected here in order to clearly illustrate 
the various types of solutions labelled in Fig. |3(b)| At lower 
temperatures the form of the velocity solutions is similar. 

At the critical radius Vc either the flow velocity equals 
the sound speed, v = Cs, or dv/ds — 0. There are several dis- 
tinct velocity solutions labelled in Fig. |3(b)| For very small 
subsonic initial Mach numbers (curve A) the flow remains 
subsonic all the way to the star, and for large supersonic ini- 
tial Mach numbers, it remains supersonic from the disc to 
the star (curve B). There is also a range of initial Mach num- 
bers (curves C) where the flow will not reach the star, and 
one value of M that results in a transonic solution - where 
the flow leaves the disc at a subsonic speed and accelerates 
hitting the star at a supersonic speed (curve E). Observa- 
tions of the widths of line profiles suggest that the accret- 
ing material reaches the stellar surface at s everal hundred 
kms~ ^, certainly at supersonic speeds (e.g. lEdwards et alJ 
I1994I) . The velocity profiles indicate that it is possible to 
have accretion flows which leave the disc at a low subsonic 
velocity but which arrive at the star with a large supersonic 
velocity. In Fig. |3(b)| the transonic solution arrives at the 
star with a Mach number of 6.54, which at a temperature 
of fO^K corresponds to an in-fall velocity of 243kms~^. For 
a realistic accretion temperature of f 0'*K the in-fall velocity 
is 259kms"^ 

Models of funnel flows h ave been studied using an 
isothermal equatio n of state llLi fc WilsonI If 99Sf) and for 
a polytropic flow iKoldoba et al.l l200^i ~with the latter 
demonstrating that transonic accretion flows are only possi- 
ble for a range of start ing radii around the coro tation radius. 
MHD simulations by iRomanova et al.l i2002f) also indicate 
that accretion flows can arrive at the star with large su- 
personic velocities. However, our aim here is to determine 
whether or not the magnetic field geometry has any affect 
on accretion; in particular how the field structure affects 
the mass accretion rate, rather than to discuss the types of 
velocity solutions that we would expect to observe. In Ap- 
pendix A we discuss an efficient algorithm which allows us 
to determine both the location of the sonic point on a field 
line and the initial Mach number required to give a smooth 
transonic solution. This algorithm may be applied to accre- 
tion flows along field lines of any size, shape and inclination, 
even in the absence of analytic descriptions of the magnetic 
field and effective gravity. 

3.3 Mass accretion rates and filling factors 

For accretion to occur the effective gravity at the point 
where a field line threads the disc must point inwards to- 
wards the star. From this subset of field lines we select those 
which are able to contain the corona and support transonic 
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accretion flows. We further check to ensure that the plasma 
beta resuhing from accretion remains < 1 along their length. 
Therefore, at any particular azimuth, accretion occurs along 
the flrst field line at, or slightly within, the corotation radius; 
the field line must be able to contain the coronal plasma and 
have a sonic point along its length. In order to determine if 
a field line can support a transonic accretion fiow, the first 
step is to find the pressure and velocity structure its length, 
which will be similar to those described in §3.2. To do this 
we need to determine the initial Mach number that would 
produce an accretion flow. We can achieve this by determin- 
ing if a field line has a sonic point as discussed in Appendix 
A. 

To calculate a mass accretion rate we require the veloc- 
ity and density of each accretion fiow at the stellar surface, 
and also the surface area of the star covered in hot spots. 
For an assumed accretion fiow temperature we determine 
the initial Mach number required to generate a transonic 
velocity profile, along each field line, and determine the in- 
fall velocity from 12011 . At every point along a field line we 
know the ratio of pressure at that point, to that at the disc, 
p/Pd- For an isothermal equation of state p cx p, so we also 
know the ratio of densities p/ pd at every point along the 
accreting field line. Thus for a given disc midplane density 
Pd, we can estimate the density at the stellar surface p,. 
Throughout we assume a constant disc midplane density of 
Pd = 5.0 X 10~®(?cm~'^, a res onable value at the corotation 
radius for T Tauri stars (e.g. iBosjllQQd) . The mass accre- 
tion rate may be expressed in terms of quantities defined at 
the disc plane, with M oc pd- Therefore raising or lowering 
Pd directly increases or decreases M. We estimate the to- 
tal surface area of the star covered in accretion hot spots 
by summing the area of individual grid cells which contain 
accreting field line foot points. For each grid cell i (of area 
Ai) on the stellar surface we obtain the average in-fall ve- 



locity Vt and average density p, of material accreting into 
that cell. Most grid cells do not contain any accreting field 
line foot points and therefore do not contribute to the mass 
accretion rate. The mass accretion rate is then the sum over 
all cells i containing accreting field line foot points. 



M = Y,Kh=Y,[Av.-p.]^ 



(28) 



The mass accretion rate can be expressed equivalently as 
M = pd^dAd, where Ad is the surface area of the disc that 
contributes to accretion (which depends on the radial extent 
of accreting field lines within the disc). Using the surface 
area of grid cells within the disc which contain accreting 
field lines to estimate Ad, we obtain M values that are com- 
parable to those calculated from I28I I. Therefore it makes 
little difference which formulation for M is used. The accre- 
tion filling factor face, the fractional surface area of the star 
covered in hot spots, is then calculated from. 



/a. 



(29) 



We have calculated the total mass accretion rate that 
dipole fields can support for isothermal accretion fiow tem- 
peratures of between 10"^K and lO'^K; for values of f} from 
0°to 90°, where f) is the obliquity of the dipole (the angle be- 
tween the dipole moment and the stellar rotation axis); for 
the DF Tan parameters in Tableland for a coronal temper- 
ature of lOMK (see Fig.gJ. = 0° and /3 = 90° correspond 
to the aligned and perpendicular dipoles respectively. 

For dipolar accretion we obtain typical mass accretion 
rates of M « lO"'' - lO-^Mej/r"!. The mass accretion rate 
increases with temperature in all cases, but at low accretion 
temperatures Tacc, there is little difference in M for all val- 
ues of P (see Fig. |4(a)j . For higher Tacc values, the aligned 
dipole field can support mass accretion rates which are a 
factor of two times less than those fields with large values 
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Figure 4. (a) The change in mass accretion rate and (b) the accretion filling factor as a function of /3, for accretion to dipolo fields 
where /3 is the obliquity - the angle between the rotation and magnetic poles, for accretion flow temperatures of lOOOK {solid), 2500K 
(dotted), 5000K (short dash), 7500K (long dash) and lOOOOK (dash-dot), (c) and (d) show how the contribution to the accretion filling 
factor from accreting closed and open field changes with /B. The DF Tau parameters from Table [Tl have been used. There are no open 
accreting field lines for the Tacc = lOOO-ft' case. 



of /3. This, in part, can be attributed to the increase in the 
amount of open field lines which thread the disc, and are able 
to support accretion as /3 is increased (see Fig. |4(d)] l. As the 
dipole is tilted from /3 = 0° to 10° the mass accretion rate 
is reduced (see Fig. |4(a)j . This can be understood by the 
changing shape of the closed field lines as /3 is increased. For 
accretion along aligned dipole field lines, accreting material 
may flow along two identical paths from the disc to the star; 
that is it may accrete either onto the northern, or the south- 
ern hemisphere. Once the dipole has been tilted through a 



small angle, the path along the field onto each hemisphere 
changes, with one segment of the closed field line loop being 
shallower than before and curved towards the star, and the 
other being longer. This longer segment bulges out slightly, 
so that material flowing along such held lines follows a path 
which initially curves away from the star, before looping 
back around to the stellar surface. This creates a difference 
in initial Mach numbers necessary to create transonic accre- 
tion flows along the different fleld line segments, with the 
net result that some closed field line segments are no longer 
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Figure 5. The stellar surface with white (black) points indicating the closed (open) accreting field line foot points for accretion to a 
dipole with obliquity (a) (3 = 0°, the aligned dipole, where accretion proceeds onto two rings in opposite hemispheres; (b) fd = 15° where 
the accretion rings have been distorted and open field lines produce the small bands centred on 180° and 360° longitude; (c) (3 = 45° 
where accretion occurs predominantly along the open field lines and (d) (3 = 90° , the perpendicular dipole, where accretion occurs in bars 
ar ound the star's equat or. All are for an accretion flow temperature of lO^K. The average surface field strength matches that considered 
bv ljardine et all i200d) with yellow (blue) representing the positive (negative) magnetic pole. 



able to accrete transonically when /3 = 10° (see Fig. |4(c)[ l. 
As the dipole is tilted further from (3 = 10° to ^ 30° - 40° 
the mass accretion rate increases in all but the lowest Tacc 
cases. This is because once the dipole has been tilted far 
enough the open field lines (those that have foot points at 
latitudes closer to the magnetic axis) begin to intersect the 
disc (see Fig. |4(d)} . There are therefore more possible paths 
that material can take from the disc to the star, causing an 
increase in M (again in all but the lower Tacc cases). As /3 is 
further increased the amount of accreting closed field lines 
continues to reduce, whilst the amount of open field lines 
threading the disc reaches a maximum, and we therefore see 
a trend of falling mass accretion rates towards the largest 
values of /3. 



The accretion fiow temperature is important in deter- 
mining whether open field lines are able to support transonic 
accretion fiows. From Figs. |4(c')1 and |4(d)| it is clear that the 
contribution to accretion from the closed field is constant for 
all values of Tacc, whereas the contribution from the open 
field depends strongly on Tacc, with more open field lines 
accreting at higher accretion flow temperatures. At the low- 
est accretion flow temperature which we consider (lOOOK), 
there are no open fleld lines able to support transonic ac- 
cretion, even for the large values of /? where there are many 
such field lines passing through the disc. This can be under- 
stood as follows. For transonic accretion a sonic point must 



exist on a field line. At a sonic point v = Cs; applying this 
to 1201 1. substituting into 11911 and rearranging gives. 



1 dv , 2 
V as ^ 



:)=g.//-§--^, (30) 

from which it can be seen that there exists some critical 
radius Vc where either w = Cs or dv/ds = 0. Clearly at this 
critical radius the two terms on the RHS of 1301 must be 
equal. 



dB 



- Seff 



(31) 



B ds 

where all the terms are evaluated at re- It should be noted 
that 13111 may also be obtained by finding the maximum 
turning point of l|A5fl . consistent with our argument in Ap- 
pendix A. 

The condition for a sonic point to exist on any field line, 
open or closed, may be expressed as equation H31^ . § is a unit 
vector along the path of the field which may be written as 
s = B/_B, which we can use to rewrite 1311 as, 

cJ^=ge„-B. (32) 

From this it can be seen that the condition for a sonic point 
to exist depends on three things: first, the path that a field 
line takes through the star's gravitational potential well; sec- 
ond, how quickly the strength of the magnetic field varies 
as the flow moves along the field line and third, the accre- 
tion flow temperature which enters through the sound speed, 
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(a) (b) 

Figure 6. Surface magnctograms as used in Fig.^for (a) the LQ Hya-like and (b) the AB Dor-like magnetic fields, coloured to show 
the strength of the radial component of the field with red representing IkG and black -IkG. White points are the accreting field line 
foot points and give an indication of the location of hot spots. Hot spots span a range of latitudes and longitudes; this is in contrast for 
accretion to a dipole field where the accreting field line foot points would be at high latitudes towards the poles. The accretion filling 
factor is about 2% in both examples. 



where for isothermal accretion oc T. It is the interplay be- 
tween all three of these factors which determines if a sonic 
point exists. 

For low accretion temperatures the sound speed Cs is 
small, whilst for open field lines the geff ■ B on the RHS of 
13211 is usually larger than for closed field lines. This can be 
seen by considering the simple example of a closed dipolar 
field line and a radial open field line threading the disc at 
the same point Rd- For both field lines the effective gravity 
vector geff at Rd is the same. For the radial open field line 
the magnetic field vector B is aligned with geff, whereas 
for the closed dipolar field line, which threads the disc at a 
large angle, there is some angle 9 between it and B. Hence 
the scalar product gef / • B = gB cos 9 is larger for the open 
field line. Thus for open field lines the RHS of (I32t is larger 
compared to closed field lines, and therefore a higher ac- 
cretion flow temperature is required to create the necessary 
high value of Cs in order to balance the two sides of equation 
13211 . Of course, for reasons discussed above, dB/ds is larger 
for the radial field line, but it is not large enough to compen- 
sate for the geff ■ B term being so small. At higher values 
of the accretion flow temperature more open fleld lines are 
able to accrete, thus helping to increase the mass accretion 
rate even though the number of closed accreting fleld lines 
is significantly less compared to an aligned dipole field. 

The aligned dipole field has an accretion filling factor 
of just under 4% (see Fig. |4(bj| , with the 13 = 15° case hav- 
ing the largest filling factor due to the shape of the closed 
field lines threading the disc, which allows material to be 
channelled onto a larger area of the stellar surface. More ac- 
creting open field lines means a larger fraction of the star is 
covered in accreting field line foot points thus increasing the 
face at large accretion flow temperatures. However, when /3 
becomes larger as we tilt the dipole further onto its side, 
face decreases. As the fleld is tilted fewer closed field lines 
are available for accretion as they no longer intersect the 
disc (see Fig. |4(c)| l, leading to a decrease in the filling fac- 
tor (see Fig. |4(b)j . The accretion filling factor is smallest 
for the perpendicular dipole where material is accreted onto 
bars about the star's equator (see Fig. |^, compared to the 



accretion rings about the poles obtained with the aligned 
dipole. 

We therefore conclude that by considering accretion to 
tilted dipolar magnetic fields both the mass accretion ac- 
cretion rate and the accretion filling factor are dependent 
on the balance between the number of closed and open field 
lines threading the disc. If there are many open field lines 
threading the disc, then a higher accretion flow temperature 
is required in order for the open fleld to contribute to M. At 
high accretion flow temperatures there is at most a factor 
of two difference between the mass accretion rate that the 
aligned dipole can support compared to the tilted dipoles 
with large values of /3. It appears as though the structure 
of the magnetic fleld has only a small role to play in deter- 
mining the mass accretion rate, at least for purely dipolar 
fields, however, as is discussed in the next section, the mag- 
netic field geometry is of crucial importance in controlling 
the location and distribution of hot spots. 



4 ACCRETION TO COMPLEX MAGNETIC 
FIELDS 

4.1 Distribution of hot spots 

Brightness modulations have long been interpreted a s evi- 
dence for hot spots on CTTSs Ce.g. lBouvier et aLllQQSll . Hot 
spots are a prediction of magnetospheric accretion models 
and arise from channelled in-falling material impacting the 
star at large velocity. The distribution of accretion hot spots 
and their subsequent effects on photometric variability have 
been studied by sev eral authors, however, only with dipo- 
lar magnetic fields l|Wp od_et ajjJljQ^jMjalijjav 
119981: IStassun fc WoodI Il999l: iRomanova et al.l |2004^ . We 
find that complex magnetic field geometries have a large 
effect on the location of hot spots. The accreting field line 
foot points are shown in Fig. |S| which represent those field 
lines in Fig. which satisfy the accretion conditions dis- 
cussed in §3.3. These give an indication of how different 
field geometries control the shape, location and distribution 
of hot spots. We find a series of discrete hot spots which 
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Figure 7. Histograms showing the distribution of in-fall velocities assuming steady state accretion flows at lO^K, using the field 
extrapolations in Fig.0 corresponding to (a) the LQ Hya-like field and (b) the AB Dor-likc field using the DF Tau (solid) and CY Tau 
(dots) stellar parameters listed in Table 111 



span a range of latitudes and longitudes with typical ac- 
cretion filling factors of around 2%. This is quite diflferent 
to what we would expect for accretion to an aligned dipole 
field, where the accreting field line foot points would be at 
high latitudes, towards the poles. With the complex mag- 
netic fields presented here hot spots can be at high lat- 
itudes, but also often at low latitudes close to the star's 
equator. The exi stence of low latitude hot spots has also 
been predicted bv lvon Rekowski fc Brandenburd ll200d) who 
consider accretion along dynamo generated stellar field lines 
rather than along dipolar field lines. It is worth noting that 
the line-of-sight field components inferred from polarisation 
measurements made usins ; the Hel 5876A emission line by 
IValenti fc Johns-Krulll i2004l) , are well matched by a simple 
model of a single hot spot at different latitudes dependent on 
the particular star. Such observations already suggest that 
low latitude hot spots may be a common feature of CTTSs. 



4.2 In-fall velocities 

In Appendix A we discuss an analytic method for calculating 
the location of the critical radius. This allows us to deter- 
mine which of the field lines from Fig. Q have sonic points, 
and then find a transonic velocity solution, where material 
leaves the disc at a low subsonic speed but arrives at the 
stellar surface with a large supersonic speed. The accreting 
field geometry obtained when considering complex magnetic 
fields is such that there are many field lines of different size 
and shape which are able to support accretion flows. This 
results in a distribution of in-fall speeds, rather than a dis- 
crete in-fall speed that would be expected for accretion along 
aligned dipole field lines. Fig. Q shows the distributions of 
in-fall velocities for our two accreting field geometries, and 
for each set of stellar parameters listed in Tabled for an ac- 



cretion flow temperature of lO'^K and a coronal temperature 
of lOMK. In both cases the in-fall speeds are large enough 
to produce hot spots on the stellar surface. We obtain larger 
in-fall velocities when using the CY Tau parameters, which 
is a higher mass star, compared to the DF Tau parameters. 
Material accreting onto larger mass stars will experience a 
steeper gravitational potential than for accretion onto lower 
mass stars, and this combined with the larger corotation 
radii for the larger mass stars, naturally leads to greater 
in-fall speeds. 

The accretion flow temperature has a negligible effect 
on the average in-fall velocity, for any magnetic fleld struc- 
ture. We flnd that the average in-fall speed remains almost 
constant as the accretion flow temperature is varied, chang- 
ing by only a few kms~^. As the temperature is increased 
the sonic point of the accretion flow is further from the disc, 
which reduces the flnal Mach number by which material is 
arriving at the star. However this reduction in the flnal Mach 
number is caused almost exclusively by the increase in the 
sound speed at higher temperature, whilst the average in-fall 
velocity remains constant. 

There is a broader distribution of in-fall velocities when 
considering the CY Tau stellar parameters compared to the 
narrow distribution for DF Tau, which is more strongly 
peaked about a single in-fall velocity (see Fig.|7|l. We calcu- 
late the natural radial extent of DF Tau's corona Rs to be 
larger than the corotation radius (Rs = 2.57?* for the AB 
Dor-like fleld and 2.77?, for the LQ Hya-like fleld compared 
to Rco ~ 2.477?,) and therefore accretion proceeds along a 
mixture of open and closed fleld lines which thread the disc 
at corotation. It can be seen from Fig. |l(a)| that accretion 
occurs almost exclusively from the corotation radius with a 
small range of azimuths where there are neither open nor 
closed fleld lines threading the disc at corotation, and there- 
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Figure 8. 

iRebuU et al 

of data from 



The correlation betwee n mass accretion rate and stellar mass. Pr eviou sly published va l ues \n oints] are taken from 
200dt): iRebuU et al.l J2002l). iMohantv et al.l ^2005^. iNatta et alj <200(j) and iMuzeroUe et al.l ipOOsIl which is a colle ction 

( ^ullbring et al .1 Jl998h . IWhite fc GheJ J20nil) . IWhite fc Basriri2003h. I MuzeroUe et al .1 i200.'i) . ICalvet et all 12004) and 

iNatta et al.l l20o3) ~Using data from the COUP sample of accreting stars lGetmaiietaLl20o3r our accretion model yields mass accretion 
rates [large circles] that are consistent within the observed scatter. Open circles are values calculated from COUP stars with estimates 
of i?*, A4t, Prot, coronal temperature and a Call equivalent width measurement, whilst the filled circles are values calculated for stars 
which are regarded as strongly accreting with W(CaII)< — lA. All quantities have been calculated using the LQ Hya surface map and 
the higher coronal temperatures derived from the COUP data, (b) shows only those stars which are regarded as strongly accreting with 
the dashed lines indicating the range in mass covered by such COUP stars. 



fore the disc extends closer to the star. This effect is much 
greater for the CY Tau parameters as its corotation radius 
is beyond the natural extent of the corona {Rs = 2.77?, for 
the AB Dor-like field and 2.97?. for the LQ Hya-like field 
compared to Tico = 9.557?*). Most of the accretion occurs 
along radial field lines from corotation. However, at some az- 
imuths there are no open field lines stretching out through 
the disc at corotation, and therefore the inner disc at those 
azimuths is much closer to the star, with accretion occurring 
along the closed field lines loops which constitute the star's 
corona. A more complete model would take account of the 
torques resulting from allowing accretion to occur from well 
within the corotation radius. We do not account for this in 
our model, however this only occurs for a very small number 
of field lines, with the bulk of accretion occuring from coro- 
tation. The number of field lines accreting at lower velocity 
from within corotation is small compared to those accreting 
from 7?co, with the resulting distribution of in-fall speeds in- 
stead refiecting variations in the size and shape of field lines 
accreting from corotation. 



We have found that the structure of these complex 
magnetic fields and more importantly the stellar parame- 
ters (mass, radius, rotation period, coronal temperature) are 
critical in controlling the distribution of in-fall velocities. 
Consequently this will affect mass accretion rates, which we 
discuss below, and also spectral line profiles, which will be 
the subject of future work. 



4.3 Mass accretion rates and stellar mass 

To investigate the correlation between mass accretion rate 
and stellar mass we require a large sample of accreting stars 
with estimates of stellar mass, radius, rotation period and 
coronal temperature. The X-ray emission from young stars 
in the Orion Nebula has been studied by COUP (The Chan- 
dra Orion Ultradeep Project), which provides a vast amount 
of data on accreting stars with estimates of all the param- 
eters required by our mo del to calculate mass accretion 
rates. ICetman et al.l (|200^ provides an overview of the ob- 
servations, and the COUP dataset, which is available from 
ftp://ftp.astro.psu.edu/pub/gkosta/COUP-PU BLIc71 

iRebuU et alj (l200Cf) first noted that the apparent in- 
crease in mass accretion rate with stellar mass and that the 
lack of low mass stars with hig h M was a real effec t for 
stars in the Orion flanking fields. Iwhite fc Ghed i200ll) also 
noted an apparent M — A7* correlation for sta rs in Taurus- 
Aurig a, with a large scatter in M values, with lRebuU et all 
i2002f) reporting that the correlation also existed for stars 
in NGC 2264. The correlation was then found to extend 
across several orders of magnitude in mass with the detec- 
tion of accretion in low ma ss T Tauri stars and bro wn dwarfs 
JWhite fc B asriF2003^ with lMuzeroUe et alj (1200311 being the 
flrst to suggest tha t M oc M^. Subsequent observations by 
ICalvet et alj (|200^ indicated that this correlation extended 
to the higher mass, intermediate mass T Tauri stars with 
several authors addi ng data at lower masses from various 
star forming regions (INatta et al .120041: iMohantv et alfeo'osi : 
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iMuzerolle et al.ll2005^ ■ with lNatta et al.l J2006l) adding data 
from p-Ophiuchus. There can be as much as three orders 
of magnitude scatter in the measured mass accretion rate 
at any particular steUar mass. It should also be noted that 
mass accretion rate measurements for stars in the Trapez- 
ium cluster are not con sistent with the M oc correlation. 
iRobberto et al.l l|20o3) report M values for the Trapezium 
cluster which are significantly lower than those obtained for 
stars in Taurus and the Orion flanking fields and suggest 
that the probable cause is that the discs of lower mass stars 
are being disrupted by UV radiation from the Trapezium 
OB stars, causing a large drop in mass accretion rates. Also, 
ICalvet et alj (l2004fl point out that there is a strong bias 
against the detection of intermediate mass T Tauri stars 
(M, = 1.5 — 4M0) with lower mass accretion rates, there- 
fore the expone nt of the M — Mt corre lation may be less 
that 2. Further, IClarke fc Pringld J200d) demonstrate that 
currently available data are limited by selection effects at 
high M values (whereby accretion rates cannot be deter- 
mined when accretion luminosity is greater than stellar lu- 
minosity) and also by a lower bound defined by the upper 
limits of non-detections. Their work therefore suggests that 
the steep correlation between M and M* is a natural conse- 
quence of detection/selection limitations, and that the true 
Af — Af« correlation may be different. 

Using the complex magnetic fields discussed above and 
assuming an accretion flow temperature of lO'^K, we have 
calculated mass accretion rates and accretion filling fac- 
tors for the COUP sample of stars which have estimates of 
M«, Rt, Prot, coronal temperature and measurements of the 
Call 8542 A line, using the lower coronal temperature, and, 
for those stars with spectra fitted with a two-temperature 
model, using the higher coronal temperature as well. We 
have looked for a correlation between the mass accretion 
rate and stellar mass of the form Af oc Af", where a is a 
constant. In Fig. |8(a)| we have plotted our calculated mass 
accretion rates from the COUP stars as a function of stel- 
lar mass, over-plotted on published values, for the LQ Hya- 
like magnetic field us ing the higher coronal temperatures. 
IPreibisch et al.l ll2005l) describe how the equivalent width of 
the Call 8542A line is used as an indicator of accretion for 
stars in the COUP dataset. They follow the classification 
discussed bv lFlaccomio etlil] (|20^) , who assume that stars 
are strongly accreting if the Call line is seen in emission 
with W(CaII)< -lA. Stars with W(CaII)>lA are assumed 
to be either weak or non-accretors. However, it should be 
noted that many of the stars in the COUP dataset have 
W(CaII)=OA, and as su ch cannot be identifie d as either ac- 
creting or non-accreting. Istassun et al.l i2006l) discuss, with 
particular reference to the COUP, how some stars can show 
clear accretion signatures in Ha but without showing evi- 
dence for accretion in Call. Therefore the sample of COUP 
stars considered in Fig. |8(a)'| may not be restricted to actively 
accreting CTTSs and could also include other non-accreting 
young stars, such as weak line T Tauri stars, where the disc 
is rarefied or non-existent, and stars which are surrounded 
by discs but are not actively accreting at this time. How- 
ever, we consider all of the available stellar parameters from 
the COUP dataset in order to demonstrate that our model 
produces a similar amount of scatter in calculated M values. 
In Fig. |8(b)| we have only plotted those COUP stars which 



are regarded as strongly accreting based on the equivalent 
width of the Call 8542 A line. 

By fitting a line to the filled circles in Fig. |8(b)| we 
find a correlation of the form M oc M^'^ for COUP stars 
which are regarded as strong accretors, using the LQ Hya- 
like magnetic field. For the AB Dor- like field we find M oc 
Ml-^. Therefore our simple steady state isothermal accre- 
tion model produces an increase in mass accretion rate with 
stellar mass, and predicts Af values which are consistent 
within the observed scatter, but it underestimates the ex- 
ponent of 2 obtained from published values. We find similar 
results when considering both the lower and higher coronal 
temperatures. However, the strongly accreting COUP stars 
only include stellar masses of A/« « 0.1 — 2.6Af0, and as 
such do not provide a large enough range in mass to test 
if our accretion model can reproduce the observed correla- 
tion. Therefore if we only consider the observational data in 
the restricted range of mass provided by the COUP dataset 
(indicated by the vertical dashed lines in Fig. |8(b)| l then 
the exponent of the observed correlation is less than 2, with 
M oc Af^'^, and if higher mass stars do exist with lower 
mass accretion rates, co uld even be less than 1.4. Therefore 
our simple model, which ITardine et al.l feOOd) have success- 
fully used to explain the observed correlation between X-ray 
emission measure and stellar mass, comes close to reproduc- 
ing the stellar mass - mass accretion rate correlation. It is 
also worth noting that our model produces a correlation in 
agreement IClarke fc Pringld J200d) . 

The accretion filling factors are typically around 2.5% 
for the sample of COUP stars (see Fig.|^, although it varies 
from less than 2% to greater than 4%. There is a slight trend 
for higher mass stars to have smaller filling factors, that is 
there is a slight trend for stars with larger corotation radii 
to have a smaller filling factors despite their larger mass ac- 
cretion rates. For stars with smaller corotation radii (lower 
mass) there are many field lines threading the disc and able 
to support accretion, and therefore a large fraction of the 
stellar surface is covered in accreting field lines foot points, 
whereas for stars with larger corotation radii (higher mass) 
there are less accreting field lines and therefore smaller fill- 
ing factors. However the higher mass stars are accreting at 
a larger velocity therefore producing larger mass accretion 
rates. The actual accretion filling factor depends on the mag- 
netic field structure. For the AB Dor-like magnetic field the 
accretion filling factors are smaller at around 1.8%. 



5 SUMMARY 

By considering accretion to both dipolar and complex mag- 
netic fields we have constructed a steady state isothermal 
accretion model were material leaves the disc at low sub- 
sonic speeds, but arrives at the star at large supersonic 
speeds. We find in-fall velocities of a few hundred kms~^ 
are possible which is consistent with measurements of the 
redshifted absorption compone nts of inverse P-Cygni pro- 
files (e.g. lEdwards et al.|[l993) . We found that for accre- 
tion along aligned and perpendicular closed dipole field lines 
that there was little difference in in-fall speeds between the 
two cases. However, for tilted dipole fields in general (rather 
than just considering a single closed field line loop) there are 
many open and closed field lines threading the disc provid- 
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Figure 9. The accretion filling factors are small, and typically 
around 2.5% for the LQ Hya-like magnetic field. Open and filled 
circles represent the same set of stars in Fig. |S] 

ing difTerent paths that material could flow along from the 
disc to the star. The path that a particular field line takes 
through the star's gravitational potential combined with the 
strength of magnetic field and how that changes along the 
field line path, are all important in determining whether or 
not a field line can support transonic accretion. At low ac- 
cretion fiow temperatures open field lines are typically not 
able to support transonic accretion, but at higher accretion 
fiow temperatures they can accrete transonically and conse- 
quently we see an increase in the mass accretion rate. We 
only find a factor of two to three difference in M values for 
accretion to tilted dipolar magnetic fields suggesting that 
the geometry of the field itself is not as significant as the 
stellar parameters (which control the position of corotation) 
in controlling the mass accretion rate. 

The magnetic field geometry is crucial in controlling the 
location and distribution of hot spots on the stellar surface. 
For accretion to complex magnetic fields we find that hot 
spots can span a range of latitudes and longitudes, and are 
often at low latitudes towards the star's equator. We find 
that the accretion filling factors (the fractional surface area 
of a star covered in hot spots) are small and typically around 
2.5%, but they can vary from less than 1% to over 4% and 
rarely to larger values. This is consistent with observations 
which suggest small accretion filling factors and the infer- 
ence of hot spots at various latitudes dValenti fc Johns-Krulll 
|2003). 

For accretion with complex magnetic fields there is a 
distribution of in-fall speeds, which arises from variations 
in the size and shape of accreting field lines. The resulting 
effect that this will have on line profiles will be addressed in 
future work. In our model most of the accretion occurs from 
the corotation radius, but at some azimuths the disc extends 
closer to the star meaning that a small fraction of field lines 
are accreting material at lower velocity. Lower mass stars, 
with their lower surface gravities, typically have larger coro- 



nae which would extend out to corotation llJardine et alj 
l200(i) . Therefore the lower mass stars are accreting along 
a mixture of open and closed field lines from the corota- 
tion radius. In contrast higher mass stars, with their higher 
surface gravities, have small compact coronae, and so the 
star is actively accreting along mainly open field lines from 
the corotation radius. However, such open field lines do not 
thread the disc at all azimuths, with some accretion instead 
occuring along field lines which are much closer to the stel- 
lar surface. This gives rise to a small peak at low velocity 
in the in-fall velocity distribution. However, this represents 
only a small fraction of all accreting field lines which do not 
contribute significantly to the resulting mass accretion rate. 

Finally we applied our accretion model to stars from 
the COUP dataset which have estimates of the stellar pa- 
rameters and measurements of the equivalent width of the 
Call 8542A line, which is seen in emission for accreting 
stars. For the complex magnetic fields we calculated mass 
accretion rates and accretion filling factors as a function 
of stellar mass. The ob served stellar mass - accretion rate 
correlatio n is M (x f'MuzcroU e et alj|2003l: ICalvet et alJ 
12004: Mohantv et al..,2005 : Muzcrolle et al.ll2005r) . however, 
this may be strongly infiuen ced by detection/selection ef- 
fects JClarke fc Pringlell2006ll . By only considering observa- 
tional data across the range in mass provided by the COUP 
sample of accreting stars, the observed correlation becomes 
A4 oc Af^ *. Our steady state isothermal model gives an ex- 
ponent of 1.1 for the LQ Hya-like magnetic field and 1.2 
for the AB Dor-like field with similar results for both the 
high and low coronal temperatures, with the caveat that the 
observed correlation may be less than 1.4 due to a strong 
bias against the det ection of higher m ass stars with lower 
mass accretion rates dCalvet et aljl20o3l . It may be the case 
that an exponent of 1.2 compared to the observed 1.4, repre- 
sents the best value that can be achieved w i th a s teady state 
isothermal accretion model. Ijardine et alJ i2006t) have used 
this model to reproduce the obs erved increase in X-ra y emis- 
sion measure with stellar mass dPreibisch et alJl2005^ ■ How- 
ever, they find that when using the complex magnetic fields 
presented here (extrapolated from surface magnetograms of 
the young main-sequence stars AB Dor and LQ Hya) they 
slightly under estimate the emission measure-mass correla- 
tion. When they use dipolar magnetic fields, which repre- 
sent the most extended stellar field, they slightly over es- 
timated the correlation. This suggests that T Tauri stars 
have magnetic fields which are more extended than those 
of young main sequence stars, but are more compact than 
purely dipolar fields. 

Also, our model has not been tested across a large 
enough range in mass to make the comparison with the ob- 
served correlation of M oc . However, it already compares 
well w ith the alternative suggestion of IClarke fc Pringlg 
i2006f) that the M — M* correlation is not as steep with 
M oc Mi'^^ . In future we will extend our model to consider 
accretion from the lowest mass brown dwarfs, up to inter- 
mediate mass T Tauri stars. This will require estimates of, 
in particular, rotation periods and coronal temperatures. 

Our accretion model reproduces a similar amount of 
scatter in calculated M values compared with observations. 
This can be attributed to different sets stellar parameters 
changing the structure of the accreting field. The stellar pa- 
rameters control the location of the corotation radius, whilst 
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it is the position of corotation relative to the natural coro- 
nal extent, which determines whether or not accretion oc- 
curs predominantly along the open field. However, for the 
moment we are restricted to using surface magnetograms of 
young main sequence staxs which may not necessarily rep- 
resent the true fields of T Tauri stars. It could be possible 
that variations in magnetic field geometry from star to star 
are responsible for the observed large scatter in the mass ac- 
cretion rate at any particular stellar mass. If there is a large 
difference in the structure of the accreting field from star to 
star this could lead to large variations in mass accretion rates 
onto different stars. The spectropolarimeter at the Canada 
France Hawaii telescope, ESPaDOnS: Eche lle SpectroPolari- 
metric Devic e for the Observation of Stars dPetit et all200 j : 
iDonatilliool) . will in the near future allow the reconstruc- 
tion of the magnetic field topology of CTTSs from Zeeman- 
Doppler imaging. This will allow a number of open questions 
to be addressed about the nature of the magnetic fields of T 
Tauri stars and allow us to test if our accretion model can 
indeed reproduce the stellar mass - accretion rate correla- 
tion. 
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APPENDIX A: ALGORITHM FOR FINDING A 
SONIC POINT AND THE INITIAL MACH 
NUMBER REQUIRED FOR A SMOOTH 
TRANSONIC SOLUTION 

In order to determine the initial Mach number which results 
in a transonic accretion flow we need to understand exactly 
what the different solutions in Fig. |3(b)| represent. These 
curves are obtained from the pressure solutions plotted in 
Fig. |3(a)1 using equation l|20|l . These pressure profiles track 
how the roots of 11911 change as the flow moves along a field 
line. The pressure function H19^ can be written in the form, 
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where a and b axe constants at any fixed point along a field 
line with, 
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For a Mach number prescribed to a flow leaving the 
disc, the pressure function has two roots, which consequently 
yield two velocities. One of the roots represents the true 
physical solution, and the other is a mathematical solution 
of no physical significance. The branches in Fig.|3(b)1 occur 
in pairs, with an individual branch tracking how one of the 
pressure roots changes as the flow travels along the field line. 
For example, suppose we had chosen an initial (subsonic) 
Mach number which resulted in a purely subsonic flow from 
the disc to the star (curve A in Fig. |3(b)] l. Then one of the 
roots, the one which gives a subsonic solution everywhere 
along the flow, is the true physical one. The second root of 
(Al) produces a supersonic solution (curve B), which is of 
mathematical interest, but has no physical meaning. Con- 
versely, had we chosen the initial Mach number which cor- 
responded to the same supersonic branch, then this purely 
supersonic branch is the one with physical meaning, with 
the other root then generating the purely subsonic branch, 
the mathematical artifact. Thus the purely subsonic and su- 
personic solutions exist provided our pressure function l|Alfl 
has two real roots when evaluated at each point along a field 
line. 
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Figure Al. The location of the critical radius rc as a function 
of a.ccretion flow temperature Tacc for starting radii i?d about 
the corotation radius (upper line, = 7.0i?«, middle line R|j = 
6.0-R* ~ Rco and lower line = 5. OR*). At lower temperatures 
rc is closer to the inner edge of the disc, and the difli'erence between 
the aligned and perpendicular dipolcs is larger. Variation of the 
starting radius has little effect for the perpendicular dipolc, and 
is of minor significance for the aligned dipole. 



At small values of p/pd the second {p/pd)~^ term in 
lAlll dominates with f{p/pd) — > oo as p/pd — > 0. At larger 
values of p/pd the logarithmic term dominates. Thus IIAH 
has a minimum which occurs when, 
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For this value of p/pd the pressure function reduces to. 
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which we will refer to as the critical function, fc- In order for 
two real roots to exist at some point r along a field line then 
fc{r) < must hold. Thus for purely subsonic/supersonic 
solutions to exist then /c < at each point along a field 
line, with /c reaching its highest value at the critical radius 
(where the two roots of the pressure function, and there- 
fore the velocity roots, are closest together). Another pair 
of solutions are the discontinuous ones labelled C and D in 
Fig. |3(b)| For these solutions /c < for r\ < r ^ Rd and 
Rt ^ r < r2- However for the domain r2 < r < ri, fc > 0, 
and therefore the pressure function has no real roots, and 
we cannot satisfy 11911 : hence there are no velocity solutions 
for these values of r. For r = ri and r = r-2, fc ~ 0, and at 
these points the pressure function has a single repeated root 
that coincides with the pressure function minimum. 

The final pair of solutions are the transonic ones labelled 
E and F in Fig. |3(b)| Here /c < for all r, except at the 
sonic point where r = Vc and fc ~ 0. Hence, for a transonic 
solution the maximum value of the critical function fc is zero 
at the critical radius Vc, but less than zero at all other points 
along the field line. This then gives us a robust method for 



Mass Accretion onto T Tauri Stars 17 



finding botli the critical radius and the initial Mach number 
required for a smooth transonic solution (provided such a 
solution exists). 

At the critical radius the minimum of our pressure func- 
tion is at its highest value; that is /c has a maximum turning 
point. Hence to determine if a transonic solution exists we se- 
lect any initial Mach number Ml and calculate how fc varies 
as we move along a field line from the disc to the star. If a 
sonic point exists on that field line, fc will have a distinct 
maximum at some point. The advantage of this method is 
that it is possible to recover the initial Mach number that 
would result in a smooth transonic flow, simply by varying 
Ml until /c = at r-c; that is we change M until the maxi- 
mum value of fc occurs at zero. This algorithm is an efflcient 
method for quickly determining both the critical radius and 
the initial Mach number which results in a transonic accre- 
tion flow. Of course, not all field lines have a sonic point (see 
§3.3) as the critical radius may either be interior to the star 
or exterior to the starting radius. This algorithm may be 
applied to accretion flows along field lines of any size, shape 
and inclination, even in the absence of analytic descriptions 
of the magnetic field and effective gravity. 

The location of the critical radius (or the sonic point 
for a transonic solution) determines the velocity with which 
material impacts the stellar photosphere. We have found the 
critical radius from the maximum value of fc for both the 
perpendicular and aligned dipoles for a CTTS with a mass 
of 0.5Mq, radius 2R0 and a rotation period of 7 days. For 
our two dipole cases the B/Bd and ge// • § terms that con- 
tribute to 1A5I I. have analytic forms given respectively by 
mill and either 1211 1 for the perpendicular dipole, or 12511 
for the aligned dipole; in both cases § = B/B. The loca- 
tion of rc therefore changes as we vary both the starting ra- 
dius Rd and also the temperature (which enters through the 
sound speed). It should be noted, however, that for certain 
choices of parameters sometimes fc has no distinct maxi- 
mum turning point, indicating that either the critical radius 
is interior to the star (r^ ^ R*), or beyond the starting 
radius {rc > Rd)- In these cases flows that leave the disc 
at a subsonic speed remain subsonic all the way to the star, 
and likewise supersonic flows remain supersonic at all points 
along a fleld line. 

Fig. lAll shows how the critical radius varies for a range 
of temperatures and starting radii for both the perpendic- 
ular and aligned dipoles. In both cases, the critical radius 
moves towards the inner edge of the disc with decreasing 
temperature. In other words as the accretion flow temper- 
ature decreases the sonic point moves along the field line 
away from the star, closer to where the flow leaves the disc. 
It is straight forward to explain why this should happen by 
considering a transonic flow which would leave the disc at a 
subsonic speed. As the temperature drops the sound speed 
decreases (cs oc Vt); therefore a flow leaving the disc does 
not have to accelerate for as long to reach its own sound 
speed and becomes supersonic sooner. Thus Vc moves closer 
to Rd with decreasing temperature; and for a high enough 
temperature the critical radius is interior to the star, and 
conversely for a low enough temperature the critical radius 
is beyond the starting radius. However, at least for this par- 
ticular set of stellar parameters, the actual change in critical 
radius location with temperature is small. 



We also found that as the critical radius moved towards 
the inner edge of the disc, the Mach number with which the 
flow arrived at the star increased. This would be expected 
however as the critical radius moves away from the star due 
to the sound speed decrease, which would naturally increase 
the flow Mach number at the star (for transonic and purely 
supersonic flows). 

It can be seen from Fig. lAll that for the perpendicular 
dipole varying the starting radius of the flow, Rd, has little 
effect on the location of the critical radius. Therefore chang- 
ing the size of equatorial field lines has a negligible effect 
on the velocity with which accretion flows (along the field) 
reach the star. For the aligned dipole changing where the 
field lines thread the disc does have an effect on the critical 
radius location and therefore on the final velocity with which 
material impacts the star. However, due to the overall small 
change in the critical radius location evident in Fig. lAll the 
in-fall velocity only varies by around 6 kms~^. It therefore 
appears that the physical size of field lines, in a dipole ac- 
cretion model, is of little importance for poloidal field lines 
(aligned north-south in the star's meridional plane), and is 
negligible for equatorial fleld lines (aligned east-west in the 
star's equatorial plane). 

The effect of the field orientation is also of little impor- 
tance for closed dipolar field lines. If we have a closed field 
line in the equatorial plane, with a maximum radial extent 
Rd, then changing the inclination of that field line so that it 
now lies in the meridional plane, has an effect on the criti- 
cal radius location (see Fig. lAH . The only major difference 
between the two cases enters through the effective gravity, 
which only has a radial component in the equatorial plane, 
but both r and 9 components in the meridional plane. How- 
ever, the difference in in-fall velocities is again very small, 
at only a few kms~^. This suggests that for accretion along 
closed dipolar fleld lines, where material is leaving the disc 
at a fixed radius, the field geometry has little effect on in-fall 
velocities; however, as discussed in §4.2, this result does not 
hold when we consider more complicated multipolar fields. 
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